Mon. Not. R. Astron. Soc. 000, 000-000 (201 1) Printed 8 July 201 1 (MN WT^ style file v2.2) 

Evidence for mild deviation from power-law distribution of electrons 
in relativistic shocks: GRB 090902B 

R. Barniol Duran^'^* and P. Kumar^* 

^Department of Physics, University of Texas at Austin, Austin, TX 78712, USA 
^ Department of Astronomy, University of Texas at Austin, Austin, TX 78712, USA 

Accepted; Received; in original form 201 1 March 25 

o 

CN ■ ABSTRACT 



Many previous studies have determined that the long lasting emission at X-ray, optical and 
radio wavelengths from gamma-ray bursts (GRBs), called the afterglow, is likely produced 
by the external forward shock model. In this model, the GRB jet interacts with the circum- 
stellar medium and drives a shock that heats the medium, which radiates via synchrotron 
emission. In this work, we carried out a detailed analysis of the late time afterglow data of 
Py ' GRB 090902B using a very careful accounting of the Inverse Compton losses. We find that 

nn I in the context of the external forward shock model, the only viable option to explain the 

X-ray and optical data of GRB 090920B is to have the electron energy distribution deviate 
^ I from a power-law shape and exhibit some sUght curvature immediately downstream of the 

shock front (we explored other models that rely on a single power-law assumption, but they 

all fail to explain the observations). We find the fraction of the energy of shocked plasma in 

magnetic field to be ^ 10^^ using late time afterglow data, which is consistent with the value 

j^ [ obtained using early gamma-ray data. Studies like the present one might be able to provide a 

link between GRB afterglow modeling and numerical simulations of particle acceleration in 
collisionless shocks. We also provide detailed calculations for the early ( < lO'^ s) high energy 
h>. ' (> 100 MeV) emission and confirm that it is consistent with origin in the external forward 

^sO ! shock. We investigated the possibility that the ~ 10 keV excess observed in the spectrum 

CO . during the prompt phase also has its origin in the external shock and found the answer to be 

negative. 
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1 INTRODUCTION 



$H The external forward shock model (see, e.g., Rees & Meszaros 1992; Meszaros & Rees 1993, 1997; Paczynski & Rhoads 1993; Wijers, Rees 
. . . & Meszaros 1997; Sari, Piran & Narayan 1998; Dermer & Mitman 1999) has proven to be a very useful concept in the study of GRBs (for a 
review of this phenomenon see, e.g., Piran 2004). The relativistic GRB jet or outflow interacts with the surrounding medium of the progenitor 
star (circumstellar medium or CSM) and drives a forward shock that accelerates the particles in the CSM, which radiate via synchrotron and 
Inverse Compton mechanisms. By modeling GRB afterglows one can learn more about the CSM medium properties (its density), general 
properties of the outflow (total kinetic energy in the shocked medium) and some details of the shock process (fractions of total energy in 
the shocked fluid imparted to electrons and magnetic fields). Particles are likely accelerated in collisionless shocks by the Fermi process and 
the resultant electron energy distribution is expected to be a decaying power-law in electron energy with index p (see, e.g., Krymskii 1977, 
Axford, Leer & Skadron 1978; Bell 1978; Blandford & Ostriker 1978; Blandford & Eichler 1987; Gallant & Achterberg 1999; Achterberg et 
al. 2001; Sironi & Spitkovsky 201 1). Various studies have attempted to calculate the value of p and to test for its universality among all bursts 
(see, e.g., Shen, Kumar & Robinson 2006). Even though in the majority of previous studies a single value of p is assumed for a particular 
burst, there is the possibility that the distribution function might deviate from a single power-law (see, e.g., Li & Chevalier 2001). 

In this paper we study the late ( > 0.5 d) afterglow of GRB 090902B in the context of the external forward shock model. In Section 
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2, we study the afterglow data for this GRB. In Section 3, we present alternatives to the scenario proposed in Section 2. In particular, we 
devote most of this section on a very detailed calculation of Inverse Compton losses, which might mitigate some of the issues presented 
in Section 2. In Section 4, we show that the only viable solution seems to be that the electron energy distribution is a little steeper for the 
X-ray radiating electrons than for the optical electrons in the external forward shock, that is, there is some curvature in the electron energy 
distribution spectrum. In Section 5, we use the results of the previous section to study the early (~ 50 s) high-energy gamma-ray data of 
GRB 090902B and find it consistent with the same origin as the late time afterglow emission. In Section 6, we explore the possibility that 
the additional power-law spectral component found at ~ 10 s, in addition to the prompt Band function, has also an external forward shock 
origin. We discuss our results in Section 7 and present our Conclusions in Section 8. 



2 LATE TIME AFTERGLOW DATA OF GRB 090902B 

GRB 090902B (Abdo et al. 2009) has been studied extensively for its high energy emission during the prompt phase detected by the LAT 
(Large Area Telescope) onboard of the Fermi satellite (see, e.g., Asano, Inoue & Meszaros 2010, Feng & Dai 2010, Ghisellini, Ghirlanda 
& Nava 2010, Kumar & Barniol Duran 2010, Toma, Wu & Meszaros 2010, Liu & Wang 2011, Zhang et al. 2011, Zhao, Li & Bai 2011, 
Zou, Fan & Piran 2011). In this paper we focus on the late time afterglow behavior in the radio, optical and X-ray bands (McBreen et al. 
2010, Pandey et al. 2010, Cenko et al. 201 1) with the assumption that radiation in these bands is produced by the synchrotron process in the 
external forward shock. We study the afterglow data after about a day of the explosion, because the optical data previous to this epoch seems 
to be dominated by the external reverse shock (Pandey et al. 2010). Here and in the rest of the paper, we use the convention that the observed 
specific flux at a particular energy, v, is given by /^(t) ex v^^t^" , where t is the observed time and /? {a) is the spectral (temporal) index. 
The GRB 090902B observations for t > 1 d can be summarized as follows. The optical data (2 eV), detected by UVOT (Ultraviolet Optical 
Telescope) onboard of the Swift satellite, shows f3opt = 0.82 ± 0.10 and aopt ~ 0.89 ± 0.03, whereas the X-ray data (1 keV), detected by 
XRT (X-ray Telescope) also onboard Swift, shows l3^ = 0.90 ± 0.13 and a^ = 1.36 ± 0.03 (Cenko et al. 201 1). 

We assume that electrons in the CSM are accelerated to a power-law in the external forward shock model, such that the electrons energy 
distribution is given by n(e) oc e^^, where p is the power-law index. In this scenario both the spectral index and the temporal decay index 
depend on p (Sari et al. 1998), therefore, one can relate a and /3 via the "closure relations" (see table 1 of Zhang & Meszaros 2004). The 
values of the spectral and decay indices will depend also on the region where the observed frequency falls in the synchrotron spectrum. 
The synchrotron spectrum is characterized by three frequencies: the self-absorption frequency (Va), the injection frequency (i/i), and the 
cooling frequency (fc); fi and Vc correspond to the synchrotron frequencies of electrons (just downstream of the shock front) whose energy 
correspond to the minimum energy of injected electrons and to the electrons that cool on a dynamical time, respectively, and Ua is the highest 
frequency at which the system becomes opaque to synchrotron absorption. 

Trying to explain both the X-ray and optical late time data for GRB 090902B in the context of the external forward shock model is not 
straightforward. We use the common terminology "slow cooling" ("fast cooling"), for i^i < Vc {uc < Vi), and constant density CSM (wind 
medium) for s = (s = 2), where the CSM density falls off as oc R'" and R is the distance from the center of the explosion. For this GRB 
we have the following options: 

Option 1: If the optical band, Vopt, and the X-ray band, i/^ , are in the same region of the synchrotron spectrum, that is, Ui < Vopt < Vx < Vc, 
then Px ~ Popt, which is supported by the data. However, the temporal decay indices in these two bands, which should be exactly the same, 
are very different. The X-ray decay is considerable steeper than the optical one, Aq = Ox — Oopt = 0.47 ± 0.06, and this particular issue 
is crucial for the rest of the paper. 

Option 2: If the optical band and the X-ray band lie in different parts of the spectrum, for instance, Vi < Vopt < fc < Vx, then the spectral 
indices should differ by A/3 = fix — Popt ~ 0.5, which is not supported by the data. Moreover, Aa is expected to be —1/4 (wind) or 1/4 
(constant density CSM), which is also inconsistent with the data - the X-ray decays too quickly. 

Option 3: If both optical and X-ray bands lie above the cooling frequency, that is, Vi < i^c < i^opt < i-'x, the discrepancies with the expected 
data and the observations are similar to Option 1. Since the temporal decay index is independent of type of medium there is no way to 
discriminate its type (Kumar 2000). This case is analogous to the fast cooling case where Vc < Vi < fopt < Vx- 

Option 4: Any possibility where fopt < fi is ruled out by the data. For the case of slow cooling, the spectrum would be Popt = —1/3 for 
both types of medium, which is inconsistent with the observed optical spectrum. Moreover, the optical light curve should be slowly rising 
(flat) for the constant medium (wind) case, which is inconsistent with the decaying light curve. For the case of fast cooling, if Uopt < Vc < Vi, 
then Popt = —1/3, and if Uc < fopt < fi, then Popt = 1/2 and Uopt = 1/4 - for both types of medium - which is inconsistent with the 
optical spectrum and light curve. 

We can see that Option 4 faces severe difficulties, thus we will not consider it any further. We will explore Options 1-3 in detail 
throughout the paper. 

First, we try to determine the type of medium that the blast wave is running into. For Options 1 and 2, the wind medium case can be 
ruled out. The reason is that, for Option 1, both the expected optical and X-ray decay would be too steep compared with the observations. For 
Option 2, one would expect Aa = Ox — aopt ~ —1/4, that is, the optical decay should be steeper than the X-ray one, which is the opposite 
to what it is seen. For these reasons, the only viable possibilities left are Options 1 and 2 with constant CSM, and for Option 3, optical and 
X-ray fluxes are independent of the type of CSM. 
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Pandey et al. (2010) and Cenko et al. (201 1) have also analyzed GRB 090902B and they both prefer Option 2 with a constant density 
CSM. Pandey et al. (2010) suggested a value of p = 1.8 ± 0.2, however, as the authors point out, the optical spectrum and X-ray temporal 
decay are inconsistent with the observed values. They appeal to optical extinction, which would make the optical spectrum agree with the 
observed value; nevertheless, the expected X-ray temporal decay is shallower than the observed one. Cenko et al. (201 1) findp = 2.21±0.02, 
which naturally mitigates the problem Pandey et al. (2010) have with the optical spectrum, but also gives an X-ray temporal decay that is 
shallower than the observed one. They suggest radiative losses to make the X-ray decay steeper - consistent with observations - however, as 
we will see in the next section, this would also affect the optical decay, and therefore it is not a viable option. 



3 SAVING THE EXTERNAL FORWARD SHOCK MODEL 

In this section, we explore a number of different possibilities that might modify the standard external forward shock model and help us 
reconcile the theory with the observations for GRB 090902B. In particular, we are interested in mechanisms that could potentially make the 
X-ray light curve steeper than expected in the simple external forward shock model. We consider the following possibilities: (i) Radiative 
losses in the blast wave, (ii) Temporal evolution of microphysical parameters, (iii) Temporal evolution of Compton-F parameter which would 
affect only the light curve of the observing band above Vc and (iv) Curvature in the injected electron spectrum. 

These possibilities have been discussed in the literature and applied to a number of GRBs. Radiative losses in the blast wave was 
considered by, e.g., Cohen, Piran & Sari (1998), Sari et al. (1998). Also, the possibility of having the microphysical parameters vary with 
time in the external shock was proposed by Panaitescu et al. (2006). Recently, Wang et al. (2010) considered the possibility of having 
Klein-Nishina suppression weaken with time so as to increase the Inverse Compton losses and steepen the > 100 MeV light curve of Fermi 
GRBs. A curvature or steepening of the injected electron spectrum has been considered by, e.g., Li & Chevalier (2001), Panaitescu (2001), 
Panaitescu & Kumar (2001a), Grupe et al. (2010). We now explore possibilities (i)-(iv) to find out if any of these can help us understand the 
optical and X-ray afterglow data for GRB 090902B. 

3.1 Radiative losses or temporal evolution of microphysical parameters 

For a constant CSM, the specific flux for an observed band v > Ui is. given by (see, e.g., Sari et al. 1998, Kumar 2000, Panaitescu & Kumar 
2000) 




„ 1 E+^ 1 3(p-l) 



e^" n2i 4 z/ 2 '\ivi <v <Vc 

Pul 3p-2 p . W 



where ee and eg are the fractions of energy of the shocked gas in electrons and magnetic fields, respectively, t is the time since the beginning 
of the explosion in the observer frame, Ekeaso is the isotropic kinetic energy in the shocked medium, n is the density of the CSM and Y is 
the Compton-y parameter, which is the ratio of the Inverse Compton to the synchrotron loss rates. Equation ^ is valid for p > 2; closure 
relations for p < 2 can be found in table 1 of Zhang & Meszaros (2004). 

According to Option 2, for p = 2.2, the expected optical flux (Ui < Uopt < Vc) is given by oc i-0-90;/-060^ consistent with both the 
optical decay and spectrum within l-cr and ~ 2-a, respectively, while the expected X-ray flux {Vc < Vx) is oc t~^'^^u~^'^^\ consistent with 
the observed X-ray spectrum within l-a, however, inconsistent with the observed X-ray decay by more than S-a. To be consistent, the X-ray 
light curve must be steepened by oc f ^''■^^. 

If the X-ray band is above Vc then, according to the second part of equation ([TJ, the way to steepen the light curve is by appealing to 
a decrease with time of EKE,iso, ee, es or an increase of Y. The decrease of EKE,iso, ee or es will also steepen the optical light curve, 
therefore, the increase of Compton-F is the only possibility that we consider. Let us explore these arguments in detail now. 

Radiative losses make the kinetic energy in the external forward shock decrease with time. For the X-ray band {i^c < Vx) withp = 2.2, 
the observed flux is oc -EJ^*^ j^^ and, thus, Ekejso should decrease as oc f ~"'^ to steepen the X-ray value to the observed value. However, 
since the flux in the optical band (ft < Uopt < Vc) is oc -Ekb.iso' 'he optical light curve will steepen by oc t""'^'', making the optical decay 
ctopt ~ 1.16 inconsistent with the observations by more than 5-a. Therefore, radiative losses cannot save the external forward shock model 
in Option 2. 

Appealing to a temporal evolution of e^ or e^ faces similar difficulties. The reason is that at > 0.5 d, for electrons radiating at X-ray 
and optical bands, Compton-F is y > few, and it roughly behaves as V oc ej e^ ' ; the dependence of flux on e^ - and also eg - turns 
out to be similar below and above Vc- Therefore, appealing to a temporal evolution of ee or e_g cannot steepen the X-ray light curve to the 
observed value, and at the same time leave the optical light curve decay unchangeqj. 

In Options 1 and 3, the X-ray and optical bands lie in the same spectral regime, therefore their fluxes have the exact same dependence on 
energy and microphysics parameters. Appealing to a temporal change of any of these would modify both light curves exactly the same way. 

^ Conversely, also in the context of Option 2, one could fix first the X-ray band with p = 2.5, but in this case the optical light curve would need to be made 
shallower by oc t^'^* . The same arguments as above prevent time- varying EfcE,iso or microphysical parameters to achieve this. 
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Lab frame view 



(a) (b) y 






Rest frame of incoming electron 
(a) (b) 



Figure 1. Geometry for scattering of a photon by an electron. The electron trajectory is denoted by a thick solid line, while the trajectory of the incoming 
and outgoing photon is denoted by a thick dashed line. We present the scattering viewed from the lab frame (top) and viewed from the rest frame of the 
incoming electron, that is, the electron is at rest in this frame before scattering (bottom). The diagrams before and after the collision are presented in (a) and 
(b), respectively. The un-primed quantities are in the lab frame, while the primed (') quantities are in the electron rest frame before scattering. 

For this reason, neither radiative losses nor temporal evolution of microphysical parameters can explain why the X-ray light curve decays 
faster than the optical one does. 

3.2 Temporal evolution of Compton-y 

As mentioned above, another possibility could be that the Compton-F parameter increases with time as (1 + F) oc i"'^^. Since Y only 
affects the flux for v^ < v - see equation lUll - it would only affect the X-ray band, not the optical one (when i^i < Vopt < J^c), and since the 
optical data already agrees with the expected value, then this possibility is very attractive. 

The Compton-F parameter that needs to be calculated is the one for electrons radiating at 1 keV, Yx, since it is the X-ray flux at 1 keV 
that needs to be steepened. Moreover, the temporal behavior (1 + Yx) oc i"'^^ should hold during the entire period of X-ray observations, 
which start at ~ 12.5 h and extend until ~ 15 d. In the following subsections, we search the 4-D parameter space (ee, e_g, n, EKE,iso) to 
determine if there is any part of the parameter space that satisfies the condition that (1 + Kr) oc i ' ^ . Before this, we present a very detailed 
calculation of Compton-F parameter for electrons of arbitrary Lorentz factor (LF), where we take into account the effect of Inverse Compton 
and synchrotron losses on the electron energy distribution self-consistently and include the Klein-Nishina cross section and electron recoil 
effects on Compton scatterings. 

3.2.1 Calculation of Inverse Compton loss 

A general calculation of Inverse Compton loss has been carried out by a number of authors (see, e.g., Jones 1968, Blumenthal & Gould 1970, 
see also, Blumenthal 1970). Inverse Compton loss has also been calculated in the context of GRB prompt emission and afterglows (see, e.g., 
Panaitescu & Kumar 2000; Sari & Esin 2001; and more recently Bosnjak, Daigne & Dubus 2009, Nakar, Ando & Sari 2010, Wang et al. 2010, 
Daigne, Bosnjak & Dubus 2011). The calculation of Inverse Compton loss, which allows us to calculate Compton-F, is not straightforward. 
Electrons cool via Inverse Compton scattering when they interact with synchrotron photons, and the same electron population is the one 
that emits the synchrotron photons. It is a problem of feedback, because the electron population in turn depends on the cooling the electrons 
experience. In this subsection we outline the calculation for Compton-F. We include Klein-Nishina effects and also relativistic corrections 
of the outgoing energy of the Inverse Compton scattered photons. 

In the co-moving frame of the electron, the energy of an outgoing photon after the scattering off of an electron is given by (see, e.g., 
Rybicki & Lightman 1979, hereafter RL79) 

x' 
^' " l + x'(l-cos(6'i-6"))' ^^^ 
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where the subscript "1" corresponds to the outgoing photon and the superscript (') means that the quantity is measured in the co-moving 
frame of the electron before scattering. The quantity x' (x'l) is defined as x' — hv' /me(? (a;i — hv[/me(?), where v' {v[) is the frequency 
of the incoming (outgoing) photon, and h, c and rrie are Planck's constant, the speed of light and the mass of the electron, respectively. 
The angle 9' (6'i) is that of the direction of the incoming (outgoing) photon with respect to the direction of the electron's momentum before 
scattering (see Figure[T](. 

The frequency of the incoming photon in the electron rest frame before scattering can be expressed in terms of the frequency measured 
in the lab frame, v, by using the relativistic Doppler formula 



I^7e(l + I3e COS 9) 



and the transformation between the angle of the incoming photon 9' in the electron rest frame and the lab frame, 6, is given by 
sin 6' = { — I sir 



7e(l + PeCOSt 



where /?e is the velocity of the electron, v^, divided by c, 7^ — l/-y/l — /3| is the LF of the electron. 
The energy in a bundle of scattered photons in the electron rest frame is given by 



dE' 



dQ. dv 



hv' 



daxp 

dn[ 



■hvidQidt 



(3) 



(4) 



(5) 



where I'^, is the specific intensity (in units of erg s~^ cm~^ sr'^ Hz^'^), dQ' (dQ'i) is the differential solid angle of the incoming (outgoing) 
photons, daxN/dfl'i is the differential cross section (the Klein-Nishina formula) and dt' is the duration of the event measured in the electron 
rest frame. The quantity in the square bracket is the number of photons per unit time per unit area moving within solid angle dO,' and 
frequency band dv' incident on the electron. 

The component of the momentum vector of the scattered photons along the electron velocity considered above is 



J' ,r,'j'K'(^')'^^KNhv'i /)'jn'j+' 

dpg — ail dv — ; --— cosdidilidt . 

hv' dii[ c 

The scattered photon energy in the lab frame is given by (see, e.g., RL79) 
dEs = 7e(d-Bs - yedp's) = 7e(l - /3e cos9'i)dE'^. 
Using the fact that the time interval in the lab frame is dt = 7edi', one finds 

dEs ,. a ol\dEs 

and combining this last equation with equation l[5), we arrive to 



dEs 
dt 



dQ dv 



>K'ie') 






/3e COS t 



(6) 



(7) 



(8) 



(9) 



In this last equation, both daj-cN/dQ'i and v'l are functions of v' and {9'i — 9'), and we can use equation l|2j to eliminate v'^. We thus find a 
general equation that describes the energy loss of an electron due to Inverse Compton cooling 

^= fdQ'dv'de') (dQ', ^-^°°°;f /"^,^, 

dt J J l + x'[l-cos{9',-9')] dQ', ' 

where the Klein-Nishina cross section is given by (see, e.g., RL79 eq. 7.4) 

d(TKN 



(10) 



dQ', 



16^ 



+ 



and, again, one can use equation ^ to eliminate v', , which yields 



d(j_f 



dQ', 



SfJT 

16^ 



i + K'[i-cos(e; 



+ x'[l - cos(6''i - 9')] - cos^{9'i - 9') 



(11) 



(12) 



l + a;'[l-cos(6'; -6")] 

Equation dlOb is general and contains no assumptions. To simplify the calculation we now make two assumptions. The first assumption 
allows us to simplify the expression of cos{6'i — 9') the following way. The integrand of dQ', in equation l llOt depends on 9' , and both x' 
and d(TKN /dQ', also depend on 9' - see definition of x' and equations ([3} and J12t . The cross section starts to fall-off steeply with angle 
only when \9', — 9'\> IjyfxJ; and also 9' ~ 7^^ - see equation l|4j. Thus, for parameters of interest to us, where 7e > 10'^ and x' < 10^, we 
have 9', ~ 0.1 and 9' ~ 10^'^. In this case, we can approximate cos{9', — 9') ~ cos^i. We call this the "head-on" approximation, since it 
corresponds to an incoming photon moving in the same direction as the incoming electron velocity vector before the collision as seen in the 
electron rest frame. We note, however, that we cannot set ^' = in equation ^, since that would overestimate v' by a factor of ~ 2. The 
second assumption we make is that the energy density in photons in the co-moving frame of the source where the photons are generated is 
distributed isotropically, therefore, we can write 
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dndv-^cy^il + Pe cos ( 
47r 



(13) 



where u,j is the source rest frame photon energy density. To derive the last expression we have used the fact that I^/v^ is a Lorentz invariant 
quantity and that dO! — dQ,{v/v'Y . With these two assumptions we can rewrite equation dlOt as 



— ^ = -^ / duu^ I dfi{l + PenY I dfi 
and the Klein-Nishina cross section is given by 



J / 1 — /3eMl duKN 

" l + a;'(l-Mi) rfK ' 



d(JKN 3(Tt 



d/i'i 8 [l+a;'(l-At'i) 



1 



l + a.'(l-/i;) 



7- +2: (1 -^i) -^1 



where /i = cos 6, fi'i = cos 6'i, and, as defined before, x' — h^y^v^i + I3efi)/mec'^. Combining equations l ll4t and l |15t we find 



"dT ^ le'^^'^'^^ / '^'^"'' / '^'^'•^ ^ ^''^ / '^'^'^ 



1 - Ml 



+ 



x'(1-m'i) = 



Ml (1-m'i) 



[l + x'Cl-M'i)]-* [1 + :c'(1-m'i)]^ [1 + x'(1-KP 



where we made another approximation, which is that /?e ~ 1, which is valid for the case, where 7e ^ 1 studied here. 
The n'l integral in equation l |16t can be carried out analytically. Let us define a function 

^ /„\ [^ Vi 



-"'-' 7_,[i + x(i-K)]"- 


It is straightforward to show that 


r^ .^^ ln(l + 2x) 
Gi{x) = , 


r* (t) 


1 1 


'~'n\X) 

(n — Ija: 


(l + 2a:)"-i 



for n 7^ 1. 



With this last equation, we can evaluate the jj,'i integral of the first term in equation J16t , which is 



dfi'i- 



1 / 

1 - Ml 



1 d ^ , ,, 2 (3 + 2a;' 
■jG3(x - 



3(l + 2a;')3' 

ln(l + 2x') 2(1 + 3a;') 



_^ [l+x'(l-/i'J]4 3 da; 

and also the second term, which yields 

, x'(l - fJ.'if _ x S , 

., ''Mi+^'(1-m'i)]^ " 2 da;'^'"'^''^" x'^ x'\l + 2x')^ 

To carry out the integral of the third term, we define a new function 

Jnix) = / d/i'i- 



Mi 



j_i [l + a;(l -/.;)]"' 
which can be shown to be 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



Ji(a;) = 
Jn (a;) = 



^ln(l + 2x)-2^ 



1 



1 + 



-G„-i(a;), forn / 1. 



{n-l)x[^' (l + 2a:)"-ij (n - l)x 
Using partial fractions decomposition, the fi'i integral of the third term in equation ( |16t can be written as 
1 



, ^'1 (1 - /i'l) 1 , 2 + a;' , 1 + a;' , 



and by using equation ( 122b . we find 



,-/,/ ^1 (l-Mi) 
,1 [!+-'(! -m'i)P 



2a;'+3, ,, , ^ ,, 2a;'^ + 20a;'' + 22x' + I 

i — ln(l + 2a; ) H =- — 

x'^ ^ ' x'-\l + 2x')3 



Combining equations ( |19l >, ( 1201 ) and | |24I >. the fi'i integral in equation | |16I ( can be written as 

,^, ,, a;'V2a;' + 3, ,, ^ ,, 2(22x''' + 75x"' + 99a;'' + 51a;' + 9) 

K{x ) = ; ln(l + 2a; ) ^ r- -^^ -. 

^ ^ x'* V -r ; 3x''^(l + 2a;')3 



(22) 



(23) 



(24) 



(25) 
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In addition, we make one last transformation of equation l ll6l l and that is tlie following; Since x' depends on /x = cos 0, we write it as 



r(l + ^eCos6')RiX(l + ^), (26) 



where we have defined a new variable X — hv^e/mi^c? , and we again have made use of the assumption that 7e ^ 1, so that /3e ~ 1. With 
this, equation J16t can be rewritten as 

—JT ^ T^^TCj'i / du-^ / dx'x'^K{x'). (27) 

The a;' integral can be done analytically, except for one term in the integrand, which is 2 ln(l + 2x')/x' , whose integral is a polylogarithm, 
specifically a dilogarithm or Spence's function. 

The rate of Inverse Compton energy loss can be expressed as 

'^ I 2 ^ dEs ^ [, dEs ,-„, 

where dE^/dt is given by equation {27}, and the second term is the integral of photon energy before scattering which is negligible compared 
with dEs/dt for 7e > 1. 

Finally, the Compton-Y parameter, defined as the ratio of the rate of Inverse Compton energy loss and the rate of synchrotron energy 
loss, is given by 

The only thing left to determine now is the energy density in photons, Uu- Since in the case of GRBs the entire source (jet) is moving 
relativistically towards the observer, then the photon energy density needed is the one measured in the source co-moving frame. All quantities 
considered in the following calculation will be in the source co-moving frame. 

The power emitted in An sr at frequency u by an electron in the shell (in units of erg s~^ sr^^ Hz^^) is (see, e.g., RL79) 

P..li^(jL]"\ ,30, 

where e is the electron charge, B is the magnetic field in the shell co-moving frame and i^^^ is the synchrotron frequency of an electron with 
LF 7e, which is given by 

eB-f! 



2TVmeC 



(31) 



The numerical factors in the last two expressions are different than the ones in Wijers & Galama (1999) for 2 < p < 3 by only less than ~ 10% 
and ~ 40% for equations l l30t and l l31t . respectively. P^ in equation i30t is valid for v < u-,^ , for f > v^^ we take P^ to vanish, even though 
strictly speaking it decreases exponentially, which introduces a very small error. 

We will now calculate the specific intensity, /^, in the middle of the shell (see Figure |2j. Let us assume that the column density of 
electrons (emitters) in the shell is N (number of emitters per unit area) and the shell radius and LF are R and F, respectively. Since the 
photons cannot arrive at a point in the middle of the shell from a distance larger than R/T, which is the radius of the causally connected 
region, then the specific intensity in the middle of the shell is approximately given by 

{^flsec^l for < e ^ 7r/3 and 27r/3 < e < TT 
(32) 
£^ for^/3<e<2^/3, 

where we have used a rectangular slab geometry since R/T <^ R. Therefore, we find that the energy density in photons of frequency u in 
the middle of the shell is 



"-!/-M»)^^(!^.l) 



(33) 



A similar analysis shows that the photon energy density near the inner or outer edge of the shell is Ui, = P^N/2c. Thus, the average value 
of Ui, in the shell is Uu ~ 0.7P^N/c. 

The column density of electrons can be written as a function of the electron energy distribution, which is defined as n(7e) (number of 
electrons per unit area per unit 7e). We define the specific flux in the co-moving frame of the shell, fi, = P^N, thus, Ui, ~ 0.7/^/c. The 
specific flux is given by 

/, = /d7en(7e)P., (34) 
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Figure 2. Schematic of a GRB jet. The width of the shell in the co-moving frame of the shell is R/T, where R is the distance from the center of the explosion 
and r is the LF of the source. To calculate the specific intensity in the middle of the shell one needs to take into account all photons produced by electrons 
(emitters) within the causally connected region of radius R/T. Since R/T ^ ij we assume a rectangular slab geometry. 



where n(7e) is calculated self-consistently by solving a coupled set of equations for n(7e) and radiation (see below). Using equation l l30t 
we rewrite the last equation as 



u 



djefllje 



e-^B 



2/3 



where 7^ is the LF of electrons radiating at synchrotron frequency u, that is - see equation ( I31l > -, 



2 2TVmeCl' 

eB 



(35) 



(36) 



The task now is to determine the electron energy distribution, 7i(7e), in order to calculate the synchrotron flux given by equation l |35t . 
The electron distribution function is determined by solving the following continuity equation 



where the source term, Si^e), is given by 
'0 7e < 7», 



(37) 



(38) 



A'' is the total number of electrons crossing the shock front per unit time per unit area, and 7^ is the minimum LF of electrons injected in the 
shock. 

We can solve equation ([37} approximately the following way. Let us define a cooling time for electrons, t'^^oi, in the shell co-moving 
frame, as follows fcooi(7e) = 7e/7e and a co-moving dynamical time, tco> ^s t'^„ = R/{Vc). The dynamical time in the shell co-moving 
frame is related to the one in the observer frame, tabs, as t'^^ = Ftobs/i^ + z), where z is the redshift and the column density is related to 
the shock radius, R, as nR/3, therefore, A'^ = nR/3t'co- 

If ^'cooi > i'co for a given 7e, electrons of this LF have not cooled much in the available time, therefore, equation ( I37t reads 



n(7e) ^ ico5'(7e), for i'^oo, > t'^^. 

On the other hand, when t'^aoi < ^'co for a given je, then n(7e) wouldn't change with time and thus 



(39) 
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A [7,n(7e)] ~ S(7e). (40) 

o-ye 

or 

n(7e) = / djiSi-yi) / ie, lort'^^oi <t'co- (41) 



To summarize, we solve numerically the following equations simultaneously to find Compton Yi^e)'- equation ([37} to find the elec- 
tron energy distribution, equation ([35} to find the synchrotron flux, equation Ml\ to find the energy loss rate due to Inverse Compton and 
substituting all these results into equation ( 1291 ) gives y (70). 



4 SOLUTION FOR GRB 090902B 

We use the parameters of GRB 090902B, z — 1.8 (Cucchiara et al. 2009), luminosity distance of d_L ~ 4.3 x 10^* cm and we also use 
p — 2.2 (see Section 3.1) and investigate Option 2 (ui < Uopt < Vc < i^x) in detail. For this, we use the formalism developed in the previous 
section. 

We use a numerical code that scans the 4 external forward shock parameters, viz. te, es, n and EKE,i3o, and finds the subset that 
satisfies a set of chosen constraints (see below). The parameters are varied in the following ranges: te — 10^^ — 10^"'^, es = 10^^ — 10~^, 
n — 10~* — 10 cm~^ and Ekeaso = 10^'* — 10^^'* erg. The value of the isotropic energy released in gamma-rays for GRB 090902B was 
Ej^iso = 3.63 X 10^* erg (Abdo et al. 2009), this is the reason why we chose to vary EKE,iao from being 30 times smaller to ~ 10 times 
larger than E-y^iso- The radiation efficiency of the prompt gamma-ray emission is given by r; = E-y^iso/{E-,,iso + Ekeaso)- This efficiency 
will be useful when trying to further constrain our results. 

We numerically calculate the allowed subspace of 4-D parameter space for GRB 090902B by imposing only two constraints: (1) the 
theoretically calculated X-ray flux at 1 keVat 12.5 h, the first X-ray data, must be consistent within l-cr with the observed value of 0.38± 0.10 
/ijy (Cenko et al. 201 1, Pandey et al. 2010) and (2) the injection frequency (Ui) at this same time should be below the optical band (i/i < 2 
eV) so that the optical light curve decays with time for t > 12.5 h. 

For the resulting subspace of 4-D parameter space we calculate a few quantities of interest for this GRB. (1) Since we want to address 
the possibility of the X-ray decay being faster due to a time dependence of the Compton-Y parameter as {1 + Yx) oc i"'^^ (for Option 2: 
i^i < Vopt < i^c < Vx), we calculate Y^ for electrons radiating at 1 keV at two times, f 1 = 12.5 h and i2 = 10^2, which we define as Yx,i and 
Yx,2, respectively. These times span almost the entire duration of the X-ray observations for GRB 090902B. Steepening the X-ray light curve 
to the observed value would require logjQ[(l -I- Kc,2)/(1 + "^1,1)] = 0.21. (2) We calculate the synchrotron cooling frequency at t\ to see if 
it is consistent with the orderings considered above, that is, Uopt < fc < fx- (3) We calculate the optical flux at 21 h, which we compare to 
the observed value of 12.0 ± 0.1 fiJy (Pandey et al. 2010). (4) The radio flux at 8.46 GHz is also calculated at 5.6 d, and it is compared to 
the observed value of 130 ± 34 /iJy detected by the Very Large Array (Cenko et al. 201 1). (5) We calculate the Compton-K parameter for 
electrons radiating at < 2 eV (optical), Yo, at ti and f2, to see if Yo and Yx behave differently or not. The results of our numerical calculation 
are shown in Figure [3] 

As can be seen from Figure[3] no part of the parameter space reaches the desired value of logj^o[(l + Kc,2)/(1 + ^a;,!)] = 0.21. The 
maximum value reached is logio[(l + Yx,2)/{i + Yx,i)] ~ 0.03, which means that the X-ray decay would only steepen at most to i~^'^*, 
which is inconsistent with the observed value by more than 5-a. Moreover, the optical and radio fluxes are inconsistent with the observed 
values by a factor of ~ 20 and ~ 100, respectively (Fig. [5}. More importantly, we also find that Vc at 12.5 h is actually below the optical 
band, which is inconsistent with the frequency ordering we are considering here (see Fig.|3}. All these issues rule out the possibility that the 
temporal evolution of Compton-F can explain the observed X-ray data. 

In Option 1, the optical and X-ray bands are below Uc, and thus the fluxes in these bands do not depend on Compton-Y - see equation 
([T}. For this reason the temporal evolution of Compton-Y cannot be invoked in this case to reconcile the difference between the data and the 
expectation of the external forward shock model. 

Let us now explore Option 3, where both optical and X-ray bands lie above the cooling frequency, therefore the fluxes in both bands in 
this case depend on Compton-Y. Forp = 1.5, the optical flux would decay as oc ^~(3p+io)/16^-p/2 _ ^-0.91^^-0.75 ^^^ j^^^jj^ ^^^ optical 
spectrum and optical decay index are consistent with the observed values within 1-cr. The X-ray decay then, must be steepened by f"'''^^ 
to match the observed value, therefore, we check to see if (1 + y^) ex f"'^^ is allowed by a subset of the 4-D parameter space. Moreover, 
(1 -|- yo) should evolve very slowly with time, otherwise the optical light curve would also steepen and that would be inconsistent with the 
observed optical data. 

Our numerical code is unable to handle p < 2, thus we cannot calculate the allowed parameter space forp = 1.5. However, we can run 
our code for p — 2.05 and compare our results with p = 2.2. We impose the same two observational constraints as before. For p — 2.05 we 
find that the magnitude of logjQ[(l + ya;,2)/(l + yz;,i)] remains roughly the same as it is for p = 2.2; however, logjQ[(l + yo,2)/(l + yo,i)] 
increases by a constant factor of ~ 2. This means that (l + Yo) increases with time faster than (1 + yr). We expect this behavior to continue 
when we decrease the value of p to p — 1.5 (there is no reason for an abrupt change when p falls below 2) therefore, the effect of Inverse 
Compton cooling on the electron energy distribution in this case is to steepen the optical light curve more than the X-ray light curve. This 
suggests that the temporal evolution of Compton-y in Option 3 to explain the observed data can be also ruled out. 
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Figure 3. We solve for the allowed 4-D external forward shock parameter space for constant CSM and p = 2.2 (eg, eg, n, Ekeaso) by imposing only 
two constraints: (i) The X-ray flux at ti = 12.5 h should be consistent with the observed value and (ii) the injection frequency, Ui, should be < 2 eV at the 
same time (see text). For this allowed parameter space we calculate (blue points) the cooling frequency at t\ (top left), the optical flux at 21 h (top right), the 
radio flux at 5.6 d (bottom left) and the Compton-y for electrons radiating at 2 eV, Yo (bottom right). All quantities are plotted as a function of the ratio of 
Compton-y for electrons radiating at 1 keV, Yx, at two times, t\ and t2 = lOii, with subscripts "1" and "2", respectively. We only plot the data for which 
(1 -I- Yx) increases with time. In order to steepen the X-ray light curve for GRB 090902B, with i/,; < Vopt < Vc < Vx (Option 2), to make it consistent with 
the observed behavior requires logj^Q[(l -f yc,2)/(l -I- yE,i)] =0.21, which is not found for any point in the 4-D parameter space. Also, we find i/c < 2 eV, 
inconsistent with Option 2 (the horizontal dashed line shows t'c = 2 eV - top left). Moreover, the optical and radio fluxes are inconsistent with the observed 
values (horizontal dashed lines). The bottom right panel shows that (1 + yc) <x (1 -|- Yo), and therefore whenever the X-ray light curve is steepened due to 
the increase of Yx with time, the optical light curve is also steepened by the same amount. 



4.1 What is the real solution? 

If the X-ray, optical and radio data are consistent with the external forward shock model predictions, then we should be able to find a subset 
of the 4-D parameter space for which the optical and radio fluxes agree with the observed values (the X-ray flux agrees with the observed 
value by design). Although we just found out that this subspace does not exist when we require (1 + K^) to increase with time (see Figure 
[5}, this subspace does exist when {1 -\-Yx) decreases with time (not plotted in Figure|3j. 

We take now the result of the parameter search in the previous subsection and further constrain it with the following conditions: (i) The 
external forward shock optical flux at 21 h should be consistent with the observed value within 1-cr, and (ii) the external forward shock radio 
flux at 5.6 d should be consistent with the observed value within 1-cr. This gives us a subspace of 4-D parameter space for which the X-ray, 
optical and radio fluxes as predicted by the synchrotron process in the external forward shock agree with the observed values within l-a. 
We show this subspace projected on the e_g-n plane in Figure|4]and we show Vc at 12.5 h as a function of EKE,iso/Ej^iso- We compare 
our result of es with the expectation of a magnetic field that is shock-compressed CSM field with pre-shocked value of Bo- The value of es 
downstream of the shock-front resulting from the shock compressed CSM field is es « Bi)^ / {2-Knmp(?), where jrip is the proton mass and 
nrup is the CSM mass density. As shown in Fig.|4] Bq ~ 10 /iG can explain all the afterglow radiation without the need for a strong dynamo 
amplification of shock compressed field. 

We find that i^c at 12.5 h is always > 1 keV. Note that !/c ~ 1 keV - within a factor of ~ 3 - only for EKE.iso/E^^iao < 0.06, which 
would require an extremely high efficiency of > 95% in producing the prompt gamma-rays. The overwhelming majority of the parameter 
space agrees very well with the expectations of Option 1, that is, Uc is above 1 keV at 12.5 h. However, how do we reconcile the different 
temporal decay indices of the X-ray and optical light curves? This is addressed below. 
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Figure 4. Allowed es-n plane and Uc at 12.5 h as a function of EKE,iso/E^^igo (left and right panels, respectively) for constant CSM and p = 2.2 when 
the X-ray, optical and radio fluxes predicted by the external forward shock model at 12.5 h, 21 h and 5.6 d, respectively, are consistent with the observed 
values within l-cr (Constraint 1). We further nan'ow down the allowed parameter space with the following constraints: e^ > 0.2 (Constraint 2 - see Section 
5), Ej^E,iso/J^~i.iso > 0.2 (Constraint 3) and ExE.iaol E-f.iso > 1 (Constraint 4). The points are color coded according to the applied constraints (in 
parenthesis): magenta (1), blue (1 and 2), green (1,2 and 3) and cyan (1,2 and 4). The solid black lines show the expectations for the shock-compressed 
magnetic field of seed values 1, 10, 60 and 450 /xG (from bottom to top - see text). The horizontal dashed line shows Vc = ^ keV. The injection frequency at 
12.5 h is i^i ~ 0.03 eV. Notice that these figures have many more points than the ones on Figure[3] since there are many more solutions for which {1 + Yx) 
decreases with time. For the points in these figures, we find that, at most, {1 + Yx) oc t~^-^ , however, as seen in the right panel, Vx < Vc, therefore, the X-ray 
temporal decay index is unaffected. (Please, see online version for color figure.) 

4.2 Curvature In the injected electron spectrum 

The only option left to explore is the possibility that the spectrum of injected electrons exhibits some curvature, that is, that the value of p is 
not the same for all observed bands, but that it is a function of electron energy. In this scenario, we allow p to vary and determine the best 
p values - the ones that give us the best agreement within the observed uncertainties - which are consistent with: (i) the observed optical 
spectrum and temporal decay (popt), and (ii) consistent with the observed X-ray spectrum and temporal decay (px)- 

We will consider the only viable option we have found -Pi < Vopt < Vx < Vc (Option 1) - and calculate the required values of p from 
the data. From the observed aopt and Popt, we find that popt ~ 2.2 would consistently explain the optical data to within ~ 2-a, and for the 
observed values of Ox and fix we find that px — 2.8 would consistently explain the optical data to within l-cr. 

If Ui < Vopt < Vc < I'x (Option 2), then the values ofpopt and px that best fit the data are popt = 2.2 (to within ~ 2-cr) and px ~ 2.5 
(to within 3-cr). On the other hand, if Ui < Vc < Popt < Px (Option 3), then the values of poj,t and px that best fit the data are popt ~ 1-4 (to 
within l-cr) and Pa: = 2.5 (to within 3-cr). Requiring a change in p in Option 3 from popt — lAlopx = 2.5 seems unlikely, since this change 
is very similar to the change one expects due to the cooling frequency, which for Option 3 should be below the optical band. Moreover, the 
X-ray data only agrees within 3-cr. For these two reasons, we rule out Option 3. It might seem that Option 2 provides a good solution, since 
Popt ^ Px, however, according to Figure|4l f or p = 2 . 2 (popt ~ Px) there is no parameter space where Pc < Px-lfvje repeat the calculation 
forp — 2.35, which is an average of popt andpa; in Option 2, we find that Figure|4]is basically unchanged. For this reason, we rule out the 
possibility that Option 2 with a curvature in the electron distribution function can explain the observed data. Therefore, Option 1 is the only 
viable solution for the afterglow data of GRB 090902B. 

In conclusion, we find that the external forward shock model can explain the afterglow data for GRB 090902B provided that the cooling 
frequency (Po) is larger than 1 keV at ~ 0.5 d and Pi < 2 eV. We also find that in order to explain the different temporal decay indices of the 
optical and X-ray light curves, there must be a slight curvature in the electron energy distribution function, where the spectrum of injected 
electrons steepens from ex 7^^'^ to oc 7(7^*, when 7e increases by a factor of ~ 30, corresponding to electron synchrotron frequency 
increasing from optical to ~ 1 keV. This happens effectively at a LF which corresponds to synchrotron frequency i^;,- The energy spectrum 
below (above) !/;, should be oc p"'^'^" (oc p'^'^") and the light curves should decay as oc t~^'^'^ (oc t~^'^'). In the next section we determine 
this break frequency (Pt), although, one should keep in mind that this break might not be sudden but is probably gradual. 



4.2.1 Break frequency 

Using the X-ray and optical data at a specific time, we can determine the effective break frequency necessary to explain the observations. 
At 21 h, the optical flux (2 eV) is 12 fijy and the X-ray flux (at 1 keV) is 0.2 /iJy (Pandey et al. 2010). At this time, one can show that a 
single power-law spectrum with i/""'^" reconciles these two fluxes, which means that Pt should be very close to 1 keV. We estimate Pi, at 21 
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Figure 5. 100 MeV flux as a function of ee-EKE,iso/^7,iso and Vc at 50 s as a function of Ekeaso/ E-^^iao in the left and right panels, respectively (for 
p = 2.2). The dashed Hnes indicate the observed flux at 50 s and Uc = 100 MeV, respectively. The injection frequency at 50 s is i^i ~ 0.8 keV (3 keV) for 
p = 2.2 (p = 2.5). The color coding is the same as in FigureO (Please, see online version for color figure.) 



h by assuming that the spectrum below (above) ui, is v''^'^^ {v~'^'^) and find that Vb should be 0.3 keV (this validates the calculation of the 
parameter space shown in Figured at 21 h the curvature of the spectrum is just starting to be evident in the X-ray data). One expects Ub to 
decrease with time, therefore, the evidence for the proposed curvature in the spectrum should become stronger with time. 

We use the combined spectrum presented in Pandey et al. (2010) fig. 2 at 1.9 d to determine the location of Vb at later times. The optical 
flux (2 eV) at this time is 7 /iJy and the X-ray flux (at 2.88 keV) is 0.03 /iJy. Again, assuming that the spectrum below (above) Ub is u~''^'^^ 
{v~'^'^), we can determine Vb to be ~ 80 eV. The values of Vb at 21 h and 1.9 d allow us to find that it decreases with time as ~ f "^■^^, which 
is roughly the same time dependence as Ui, which exhibits Vi oc t~^'^ . 

At the beginning of the Swift XRT observations at 12.5 h, we therefore expect Pb to be ~ 0.8 keV, and at 21 h, we expect Vb to be 
~ 0.3 keV, which is the lower bound of the XRT energy range. Therefore, there should be a very small difference in the spectrum and light 
curves between the 0.3 — 0.8 keV and 0.8 — 10 keV bands during 12.5 h to about 21 h if our interpretation that there is some curvature in 
the spectrum is correct. 



5 EARLY HIGH ENERGY DATA 

Using the parameter space we have obtained in Section 4.1 from late times radio, optical and X-ray data, we can calculate the external 
forward shock flux at high energies (FermifLAT band) at early times. We choose to calculate the flux at 50 s and at 100 MeV. If the LAT 
emission has an external forward shock origin, then this flux should be consistent with the observed value of 220 nJy (Abdo et al. 2009). On 
the other hand, the GBM band flux (~ 100 keV) is dominated by the typical prompt variable source, whose origin remains uncertain; this 
emission lasts for ~ 30 s and then it exhibits a sharp decay in the flux. The external forward shock emission at 100 keV and 50 s should be 
smaller than the observed 100 keV flux, which is ~ 0.4 mjy (Abdo et al. 2009), so that the external forward shock flux does not prevent the 
observed flux to decay very rapidly (oc t"^) as it is observed. 

The calculation of the external forward shock emission at early times is not straightforward, since we find that the spectrum of injected 
electrons has a curvature. This should be taken into account when extrapolating from about a day to a few tens of seconds. Using the estimated 
evolution of Ub as oc t~^'^ (Section 4.2.1), one finds that Ub is ~ 300 MeV at 50 s. However, the temporal evolution of Ub is highly uncertain 
and we do not know if it monotonically extends to very early times. For this reason, we use two different values of p when extrapolating 
from ~ 1 d to 50 s: p = 2.2, consistent with the p used before and p = 2.5, which is the average ofpopt and px obtained in Section 4.2. The 
results of these calculations are presented in Figure[5l 

We find a large range of values for the external forward shock flux at 100 MeV and 50 s and plot this flux against eeEKE,iso/ E-y^iso\ 
we choose this parameter because the flux above Uc is roughly proportional to eeEKE,iso (see eq.lTJ and we normalize Ekeaso to E-y^iso- 
Taking e^ > 0.2 as found for many GRB afterglows by Panaitescu & Kumar (2001b) and Ekeaso/ E-y^iso > 1 (so that 77 < 0.5), suggests 
that eeEKE.iso/ E-y^iso > 0.2, and in that case the calculated 100 MeV flux is consistent with the observed value to within ~ 1-cr (for 
p = 2.2). For p — 2.5 and eeEKE,i3o/E-,,iso > 0.2 the 100 MeV flux at 50 s is similar to the one obtained for p — 2.2, but shifted 
downwards by only a constant factor of ~ 2. This is a remarkable agreement given the fact that we have extrapolated the afterglow data from 
~ 1 d to 50 s, and from radio, optical and X-ray to 100 MeV. 

It is also important to know the location of the cooling frequency at 50 s, which would allow us to determine the spectrum and tem- 
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poral decay index of the > 100 MeV light curve. For both values of p, the cooling frequency is identical, and we plot it as a function of 
EKE,iso/E-,^iso in Figure \5\ We can see that for EKE,iso/E^^iso > 1, 100 MeV <Vc1~, 100 GeV at 50 s, therefore, one expects three 
possibilities: (i) Uc ~ 100 MeV and it will very shortly fall below the 100 MeV band, and the > 100 MeV spectrum in this case should be 
consistent with being oc v~^''^, (ii) Vc is much higher than 100 MeV and will thus remain above it until 10^ s, which marks the end of the 
LAT observations and (iii) Vc > 100 MeV and it will cross the 100 MeV band during the observations. For (i), the 100 MeV light curve 
would be /^ oc t'^-^^p-^-^° for p = 2.2 or /^ oc f-i-^^z/-^-^^ for p = 2.5. For (ii) the 100 MeV light curve would be /;. oc fO-S V"°®° 
forp = 2.2 or /^ oc t^^'^'^u'"'''^ forp = 2.5, and, lastly, for (iii) the light curve will transition from (i) to (ii). 

The observed 0.1 - 300 GeV light curve for ORB 090902B decayed as t-i-5±o.i^-i.i±o.i ^^ ^^^ ^^^^ interval 25 - 1000 s (Abdo et 
al. 2009). This light curve could be explained by scenario (i) above, for p = 2.5, within ~ 1-a. This agrees nicely with the results shown 
in Figure|4j one can see that the 100 MeV flux increases almost linearly with e^EKE,iso as expected for Uc < 100 MeV and then reaches 
a plateau. The reason for this plateau is that it corresponds to v^ Z 100 MeV at 50 s, and the flux below i/^ has already been precisely fixed 
by the X-ray flux at 12.5 h (one of the constraints), therefore, the 100 MeV flux would also be fixed precisely. Shortly after 50 s, the > 100 
MeV spectrum is consistent with being above i/c as suggested by the data. Scenario (ii) above can be ruled out, because the LAT light curve 
decreases faster than predicted and also the predicted spectrum is too shallow compared with the observed one. Finally, scenario (iii) can also 
be ruled out, since a break in the LAT light curve was not detected in the data for ORB 090902B (Abdo et al. 2009). It is worth mentioning 
that there is a small steepening of the > 100 MeV light curve by, at most, ~ ^-003 when Vc < 100 MeV, due to the increase of {I + Y) for 
electrons radiating at 100 MeV. 

The flux from the external forward shock at 100 keV, which decays as ~ t~^ would dominate the observed 100 keV light curve and 
prevent it from decaying quickly as it is observed (~ t~ ) unless the external shock contribution to the 100 keV flux is much smaller than 
the observed flux. To check this, we calculate the external forward shock flux at 100 keV and 50 s and find it to be ~ 0.01 mjy (for both 
p values), which is a factor of ~ 40 smaller than the observed value. This allows the observed prompt 100 keV flux to decay rapidly as 
observed. 



6 THE ADDITIONAL POWER-LAW COMPONENT EXTENDING TO 10 KEV 

The ORB 090902B spectrum at early times displays a Band function in the sub-MeV energy range. Time-resolved spectral analysis also 
shows a significant power-law component that appears to extend from the GeV range to the lowest energies (~ 10 keV) and it is more intense 
than the Band function both for photon energy < 50 keV and > 100 MeV (Abdo et al. 2009). The Band function is usually associated with 
the prompt emission and its origin remains uncertain. We address here the question whether the extra power-law component in the spectrum 
could have an external forward shock origin. We do this at 7 s, which is the midpoint of interval b in the analysis of Abdo et al. (2009), where 
the power-law contribution is best constrained and its spectrum is /3 = 0.94 ± 0.02. The observed flux at 10 keV and 7 s is ~ 12 mJy. 

For the external shock origin of the power-law component at 7 s, the injection frequency should be below ~ 10 keV and the cooling 
frequency should be above the highest photon energy detected at that time, which is a few GeV. Using the same procedure as in the last 
section we calculate Ui, Vc and the flux at 10 keV at 7 sq. We find that the injection frequency has a very narrow range of allowed values, 
Vi ~ 15 keV and Ui ~ 60 keV for p — 2.2 and p — 2.5, respectively. The cooling frequencies at 7 s for p — 2.2 and p — 2.5 are 
identical and are just a factor of (50/7)^ " ~ 2.7 higher than the value at 50 s (see Figure|5}, that is, 300 MeV <Uc^ 300 GeV at 7 s for 
EKE,i3o/E-,,iso > 1. However, z^c at 7 s cannot be larger than a few GeV, otherwise it will stay above 100 MeV for the entire duration of 
the LAT emission, which is inconsistent with the observed spectrum during this period. In addition, the expected spectrum between Ui and Vc 
would be /3 = 0.60 and /3 = 0.75 for p = 2.2 and p = 2.5, respectively, which is inconsistent with the observed value. Finally, the expected 
flux at 10 keV and 7 s lies in a very narrow range, and it is ~ 0.3 mJy (for both p values), which is a factor of ~ 40 smaller than the observed 
value. All these arguments suggest that the power-law detected in addition to the Band spectrum at early times is unlikely produced by the 
external shock, at least in the simplest version of this model. 



7 DISCUSSION 

In this paper we have considered the late afterglow data of GRB 090902B in the context of the external forward shock model. The optical 
data is entirely consistent with this model. However, the X-ray flux decays faster than expected for the observed X-ray spectrum. We consider 
three possibilities that could have steepened the X-ray light curve falloff. 

First, we considered radiative losses or evolving microphysical parameters that might steepen the X-ray flux. However, since the optical 
and X-ray fluxes have very similar dependence on these parameters, it is not possible to steepen the X-ray light curve and at the same time not 
steepen the optical light curve, that is, obtain a solution consistent with the observations. Second, we carried out a very detailed calculation 

■^ Note that the > 100 MeV light curve starts decaying at ~ 10 s, which might con'espond to the beginning of the deceleration phase of the external forward 
shock, however, extrapolating back to 7 s does not introduce a significant error 
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of the Compton-y parameter in hopes that its increase with time could steepen the X-ray data if the cooling frequency is below the X-ray 
band. However, we find that its effect on the steepening of the light curve is extremely small and unable to explain the observed data. 

Another modification of the standard external forward shock model we considered is the possibility that the spectrum of injected 
electrons at the shock front is not a single power-law, but that it exhibits some curvature. We find that this is the only modification that can 
reconcile the theoretical expectations with the afterglow data of GRB 090902B. 

The curvature we find in the spectrum of injected electrons might not be inconsistent with the recent particle-in-cell simulations of 
particle acceleration in relativistic collisionless shocks (Sironi & Spitkovsky 2011). Simulations need to be able to run for longer to achieve 
time-scales and energy range relevant to the GRB afterglow. Future simulations should be able to explore our claim of the presence of 
curvature in the electron energy distribution. The slight curvature, invoked phenomenologically, is tied to the physics of particle acceleration, 
which suggests in this case that the higher energy electrons - those radiating X-rays ~ are accelerated slightly less efficiently than those 
radiating optical photons. The fact that we find a "downward" curvature {p^ > Popt) is reassuring, since the opposite would have been 
contrived. 

The synchrotron frequency that corresponds to the electrons Lorentz factor at which the spectrum curves should be at ~ 0.8 keV at the 
beginning of the SwiftfXKT observations. The break frequency decreases with time as ~ i~^'^^ and it falls below the XRT band at about 21 
h (Section 4.2.1). An inspection of the 0.3-1.5 keV XRT light curve shows a slight flattening during this time in this light curve compared 
to the 1.5-10 keV one, just as predicted (Evans et al. 2007). A careful analysis of the XRT data shows that the 0.3-1.5 keV data are well fit 
by a broken power-law as predicted, but the presence of this break in the light curve cannot be claimed with high statistical significance; the 
predictions of this model however cannot be ruled out due to the sparseness of the data and the fact that the theoretical decay slopes before 
and after the break (i^"'^'^ and t^^'^^) are similar (Margutti, personal communication). Similarly, the possible break in the spectrum {v^'^'^ 
to !/~ ''■''') could not be found since the errors in the spectrum are on the order of the difference between the two different indices (before and 
after the break) that one expects. Moreover, the theoretically expected change in the spectrum and decay indices appears at the lower energy 
band of XRT, which is affected by absorption. There is the possibility, however, to use the absorption calculated using the optical band to fix 
the level of absorption in the X-ray band and look for this small change in the spectrum and temporal decay indices. 

This work assumes that the late afterglow in all bands is produced by the same population of electrons. Can we abandon this scenario 
and invoke one in which, for example, the X-ray and optical data are produced by two different sources? It seems unlikely. We discuss this 
below. 

Let us start with the assumption that the X-ray flux is originated in the external forward shock, that is, by the interaction of the GRB jet 
with the CSM. The optical data might have another origin. However, since the optical flux decays as a single power-law for a long period 
of time, its origin is also some form of an external shock. This external shock does not have to necessarily be the same one that produced 
the X-ray afterglow. It could be, for instance, that a lower LF cocoon material interacts with the CSM and drives an external shock which 
produced the optical radiation. As can be seen in Section 4.2.1, the optical and the X-ray fluxes close to 1 d fall on a single power-law 
spectrum which is roughly consistent with the observed optical and X-ray spectra. If we assume that for these two shocks te, es, and n are 
approximately equal, then EKE,iso should also be the same. Why should Ekeaso for the two different and unrelated sources be the same? 
For this reason, we think that this possibility is contrived and suggest the same origin for both X-ray and optical photons. 

Even if X-ray and optical photons are produced by the same population of electrons: Could it be that optical (and radio) photons are 
produced via synchrotron while X-ray photons originate via Inverse Compton radiation scattering of synchrotron photons? This interesting 
possibility can be ruled out. At ~ 1 d, if we extrapolate the optical flux using optical spectrum to radio band we overestimate the observed 
radio flux by a factor of ~ 700. This means that there is a break in the spectrum, which would correspond to Ui, at ~ 10~^ eV. At this 
frequency, the peak synchrotron flux is ~ 1 miy. The optical depth to Thompson scattering in the external forward shock at ~ 1 d is 
Te ~ arnR ~ 10^^, while the ratio of X-ray flux to the peak synchrotron flux at the same time is ~ lO^**. Therefore, this possibility can be 
safely ruled out since the X-ray flux is at least 10'^ times larger than the maximum expected flux for the synchrotron-self-Inverse Compton 
process. 

We have also found that when we constrain the X-ray, optical and radio fluxes at ~ 1 d one finds a large allowed range for es. When we 
further constrain e^ > 0.2, as found by Panaitescu & Kumar (2001b) and take the isotropic kinetic energy in the blast wave to be larger than 
the isotropic radiated gamma-ray energy during the prompt phase (EKE,iso > E.y.iso), then the magnetic field in the source is consistent 
with being produced via shock-compressed CSM field as shown in Figure|4]and suggested by Kumar & Barniol Duran (2009, 2010), Barniol 
Duran & Kumar (2010). The required seed field - the upstream field - before compression is < 10/^G. The constraint EKE,iso > E^,.iso 
is applied in order to avoid a radiative efficiency larger than 50%; however, even when we take Eke^xso > i?7,iso/5 - corresponding to a 
radiative efficiency of 80% - we still find that the shock compression scenario holds with just a slightly larger seed field of < 60/iG (see 
Figure |4l(. 

To obtain the region of es-n plane allowed by the GRB 090902B data we have used the X-ray, optical and radio fluxes at > 1 d, since 
the external reverse shock might dominate the optical data until ~ 0.7 d (see fig. 1 of Pandey et al. 2010). There is the possibility, however, 
that the radio flux at ~ 1 d might still be dominated by the external reverse shock, since the external reverse shock at the radio band might 
decay slower than at the optical band. This is the reason why we have obtained the es-n plane using the radio flux at 5.6 d (Cenko et al. 
201 1). We expect that at this time the contribution from the reverse shock to the radio flux is negligible. Nevertheless, if instead of using the 
radio flux at 5.6 d we use the radio flux measured at 1.3 d at the Westerbork Synthesis Radio Telescope (van der Horst et al. 2009), we find 
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Figure 6. Lower limit on the upstream circumstellar medium magnetic field, Bu, as a function of the Compton-y parameter for electrons radiating at 100 
MeV at 10^ s, 1l at (Isft panel). Yi^j^x is calculated using the method described in Section 3 for the allowed subspace of the 4-D parameter space (e^ , eg , n, 
EKE.iao) for GRB 090902B afterglow data shown in Figure|4](for p = 2.2). Bu was obtained by requiring that the Inverse Compton cooling time for these 
electrons is larger than their acceleration timescale: Method 1. We also calculate a lower limit on Bu by requiring that the electrons radiating at 100 MeV at 
10"^ s are confined to the shock front: Method 2 (Bu obtained by this method does not depend on Yi^at and has an extremely weak time dependence). The 
maximum of these two values, Bu,max, is the true minimum of Bu, which for GRB 090902B is given by Method 2 (right panel) and it is plotted as a function 



of Ylat to aid in the comparison between the two panels; B 
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that our general results are not modified: the requirements for the magnetic field are the same; however, the density increases by a factor of 
~ 10; points shown in Figure|4]are shifted to the right by a factor of ~ 10 along the diagonal bands in this case. 

The discrepancy in the value of es between our work and the work of other authors, e.g. Liu & Wang (2011), arises from the fact that 
these authors assume that the cooling frequency is between the optical and X-ray frequencies. We make no such a priori assumption (the only 
assumption we make is that X-ray, optical and radio data come from synchrotron emission in the external forward shock) and, in fact, our 
very detailed analysis shows that if the cooling frequency is between the optical and X-ray bands, it is impossible to have the X-ray temporal 
decay index match the observations. Liu & Wang (2011) find a value of es ~ 10^'^ and a density of n ~ 3 x 10^* cm~^ (for the wide 
component of their two-component jet, which dominates the flux at late times), which may give the impression that their es is ~ 10^ larger 
than our found value. However, this is not so. For the es-n region where e^ > 0-2 and Eiso/E^^iao > 1, the maximum allowed value of es 
we find for this value of n is es ~ 4 x 10"'' (see Fig. 0, which is 30 times smaller than the value reported by Liu & Wang (2011). The 
discrepancy is still large, but not at the level of a factor of 10^ and - as discussed - originates from the incorrect assumption of the location 
of the cooling frequency. Note that Cenko et al. (201 1) do not find es for GRB 090902B, but instead assume equipartition and fix their value 
of es. 



7.1 Electron acceleration and the upstream magnetic field 

Recently, the detection of 100 MeV photons up to ~ 10'' s in Fermi GRBs has been used to calculate a lower limit on the upstream CSM field 
by requiring that the Inverse Compton loss timescale is larger than the acceleration timescale of electrons radiating at 100 MeV (Li & Zhao 
2011). For the subspace of 4-D parameter space (ee, es, n, EKE,iao) consistent with the data for GRB 090902B (Section 4.1) we calculate 
Compton-Y for electrons radiating at 100 MeV at 10'^ s (using p = 2.2), and use that to determine a lower limit to the upstream field (Figure 
[6l(. We calculate another lower limit to the upstream field by requiring that electrons radiating in the LAT band are confined to the shock 
front (Piran & Nakar 2010, Bamiol Duran & Kumar 201 1); we determine this lower limit for electrons radiating at 100 MeV and 10^ s (see 
equation (3) in Barniol Duran & Kumar 201 1). The true lower limit to the upstream magnetic field is taken to be the larger of these two limits 
and the result is shown in Figure [6] We find that for GRB 090902B the electron confinement requirement gives a larger upstream field for 
much of the allowed parameter space, and an upstream field of ~ 5/iG is sufficient to confine electrons producing 100 MeV radiation and to 
avoid excessive Inverse Compton losses while electrons are traveling upstream of the shock front. 

The lower limit on the upstream CSM magnetic field we find is much smaller than the one found by Li & Zhao (201 1). This discrepancy 
arises because of two points: 1. The difference between the calculation of Inverse Compton loss - Li & Zhao (2011) used a simplified 
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calculation for Inverse Compton loss whereas we have carried out an almost exact numerical calculation as described in Section 3.2.1, and 2. 
Different numerical factors in their equation of the lower limit on the upstream fielqj. 

We note that the es-n parameter space shown in Figure |4] should be slightly revised to reflect our findings in Figure[6j results in Figure 
|4]did not include the constraints on magnetic field determined in this subsection. However, the modification to the lower limit of the CSM 
field in the es-n parameter space is very small, and it is non-existent for the case when Ekeaso/ E-i^iso > 1 and e^ > 0.2, therefore, we 
have left Figure|4]unchanged. 



8 CONCLUSIONS 

We have analyzed the late time afterglow X-ray, optical and radio data for GRB 090902B in the context of the synchrotron radiation mech- 
anism in the external forward shock model. We find that a curvature in the power-law electron energy distribution is needed in order to 
provide a good fit to the late time optical and X-ray data; radiation losses, varying microphysical parameters and an increase in Inverse 
Compton losses all fail to explain the observed data. The late time afterglow fit gives an es (fraction of total energy in the shock imparted to 
magnetic fields) consistent with shock-compressed circumstellar medium magnetic field of < 10 /iG and < 60 ^G if we take the efficiency 
for producing gamma-rays to be ~ 50% and ~ 80%, respectively. 

Particle acceleration in the external forward shock allows us to set a lower limit on the upstream circumstellar medium magnetic field. 
We find that the field strength in the unshocked medium in the vicinity of GRB 090902B must be at least 2 /iG in order to produce 100 MeV 
photons at 50 s < t < 10^ s (Barniol Duran & Kumar 201 1). 

The calculation presented in this paper represents an improvement in our previous afterglow modeling due to a more precise calculation 
of Inverse Compton losses. Here, we include Klein-Nishina effects and also relativistic corrections of the outgoing energy of the Inverse 
Compton scattered photons. We also calculate the electron energy distribution self-consistently by determining the synchrotron emission and 
using it to determine the Inverse Compton losses, which in turn modify the electron energy distribution. 

The flux calculated at 100 MeV at 50 s using the external forward shock parameters obtained from the late afterglow data is consistent 
with the FermifLAT data, confirming our previous claims (Kumar & Barniol Duran 2009, 2010). We also calculated the expected external 
forward shock at 100 keV and 50 s. At this time, the observed 100 keV light curve is undergoing a fast decay (~ t~^). We find that the 
external forward shock at 100 keV and 50 s is smaller than the observed value by a factor of ~ 40, easily allowing the observed light curve to 
decay quickly (confirming our earlier results in Kumar & Barniol Duran 2010). We speculate that for some small fraction of GRBs this steep 
decay will not be seen, instead, one will see a smooth slowly decaying (~ t^^) light curve emerge after the main prompt, variable, emission 
is over. 

The origin of the high-energy emission at < 50 s (when the prompt GBM phase is active) is still a subject of debate. It is claimed that 
GRB 090902B exhibits variability in the LAT band (Abdo et al. 2009), which would make the external forward shock origin for < 50 s 
difficult, since this model in its simplest version predicts smooth light curves. However, the fact that for this GRB the > 50 s 100 MeV light 
curve connects smoothly - as a single power-law - with the < 50 s one does argue in favor of the external forward shock origin (see, however, 
Maxham, Zhang & Zhang 2010). The delay on the detection of the first 100 MeV photons with respect of the first GBM photons from this 
GRB can also be explained by the external forward shock model (Barniol Duran & Kumar 201 1). A dedicated study of the variability during 
the prompt phase in the LAT band would certainly shed more light on whether the prompt (< 50 s) 100 MeV emission of this GRB can be 
originated from the external forward shock. 

We find that the cooling frequency at ~ 50 s is ~ 100 MeV for GRB 090902B. If this is correct, one should be able to track the cooling 
frequency as it passes through the LAT band at earlier times (i^c oc t^^' ^). A presence of this behavior in the data of this and other LAT 
detected GRBs would help confirm an external shock origin for the LAT emission. 

Finally, the prompt spectrum of GRB 090902B contains a power-law component in addition to the usual Band function. We argue 
against the external shock origin of the extra power-law component at 7 s (Section 6). 
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